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It is shown that superluminal optical signalling is possible without violating Lorentz invariance 
and causality via tunneling through photonic band gaps in inhomogeneous dielectrics of a special 
kind. 
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I. INTRODUCTION 

T-H ! 

A number of recent experiments have reported the observation of electromagnetic waves propagating with velocities 
larger than c (the velocity of light in vacuum) in dispersive media |IJ , wave guides Q , electronic circuits || and in 
tunneling The experimenters have been quick to point out that these observations do not necessarily contradict 
the special theory of relativity and causality. These claims have naturally generated a controversy in the literature. 

In the case of dispersive propagation the claim is in apparent contradiction with the pioneering work of Sommerfeld 
and Brillouin 0] who clearly showed the difference between group velocity, phase velocity and signal front velocity, 
and established the result that no physical signal can travel faster than c in dispersive media. However, it has recently 
been argued that for physical signals that are of finite duration the causality principle "cause precedes effect" is 
preserved despite superluminal motion. This is because a superluminal signal travelling backward in time can never 
arrive before the primary signal is generated, thus preventing the original user changing the transmitted signal ||. 
(f) . In the case of frustrated total internal reflection (FTIR) and tunneling, the situation is quite different. It has been 
argued that in such cases the wave number is imaginary, the phase is a constant and the concept of a signal front is 
meaningless ||. Further, it has been pointed out that if the signal is narrow-band limited, there is no distortion of 
f^*) ■ the signal envelope and its delay is the same as that of its centre of gravity Q| . Since the evanescent (exponentially 
damped) component of a wave does not oscillate with distance, it does not accumulate any phase and can therefore 
propagate through the evanescent region with zero (phase) delay. It has been argued that there is empirical evidence 
of this in, for example, symmetrical FTIR in which there is no time lag between the reflected and tunneled signals 
|^| . However, it is not quite clear how a zero phase delay necessarily implies a zero signal delay. 

One source of confusion in the literature, in our opinion, is the popular use of an analogy between the Helmholtz 
and Schrodinger equations. Since Maxwell's equations in an inhomogeneous but isotropic medium reduce to the 
Helmholtz equation for a monochromatic wave in the scalar approximation, and the Helmholtz and the non-relativistic 
Schrodinger equations are formally identical, the one-dimensional process of non-relativistic quantum mechanical tun- 
neling has been used to model the optical process of transmission through a barrier || . This is obviously unsatisfactory, 
because (a) the Schrodinger evolution used is characteristically non-relativistic whereas the optical processes in ques- 
tion are intrinsically relativistic, and (b) the Helmholtz function for the electric field is real whereas the Schrodinger 
wave function is complex. It would therefore be preferable to use a reliable and consistent quantum mechanical 
formalism for photons. 

Fortunately, such a formalism exists and is based on the classic works of Kemmer || and Harish-Chandra |jTo| . 
In this formalism, the wave function for the photon, which obeys a first-order equation similar to the Dirac equation, is 
a ten component column whose first six elements (the electric and magnetic field strengths) are real functions and the 
last four are zero, and there is a conserved four-vector current associated with energy flow (not charge flow as in the 
familiar case of charged particles with a complex wave function) whose time component is positive definite and can be 
interpreted as a probability density. The phase of such a wave function is obviously not expressible as a multiplicative 
exponential factor but is rather given in the same way as in classical electrodynamics through an additive term in 
the sinusoidal function for the fields. The signal velocity can be calculated in this formalism unambiguously from the 
energy flux vector which turns out to be proportional to the Poynting vector, as one would expect. 

It is the purpose of this paper to show, using this formalism, that Einstein causal electromagnetic signals can 
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indeed travel faster than c while tunneling through a photonic band gap provided that the dielectric in the gap is 
inhomogeneous and (practically) non-dispersive. The same result will be shown to hold for classical light. 



II. THE TUNNELING SOLUTION IN ELECTRODYNAMICS 

Let us consider the usual tunneling problem with a thin non-magnetic, practically non-absorptive material with 
a band gap around the frequency w, extending from x = to x — d and the signal incident normally on it so 
that there is no dispersion. It is essentially a two-dimensional problem (in the x — y plane) expressible in terms 
of a single component of the electric or magnetic field . We will consider the case of electric polarization with 
H x = H z = 0, E x = E y = and \i = 0, e = e(x), eo = 1. The same result will hold for magnetic polarization also. 
Then Maxwell's equations can be written in the rest frame of the dielectric material in the form 

d y E z = 0, d x E z = - c d t H y (1) 

d z H y = 0, 8 x H y = C -^-d t E z (2) 

aX-%£^o (3) 

d 2 x H y + d 2 y H y + (d x lne(x))d x H y - ^d 2 t H y = (4) 

Let us first assume that the time variation of the electric and magnetic fields is given by exp(±iu;t) , and use the 
ansatz E z (x, y) = Y{x)U(y). Then it is easy to show that 

U(y) = pe^y (5) 

where (3 and a are constants. It follows from (nl) that a — 0, and so we have 



E z = pY(x)e ±luJt (6) 
Tic/3 dY_(x) n±lult 
dx 



This shows that the magnetic field H y is completely determined by the electric field E z . It also follows from (|l|) and 
(g) that 

d -m = -^J <X )Y( X ) dx (8) 

or, 

d 2 Y(x) uj 2 r . , s 

An approximate solution to this equation (^|) is given by 

Y(x) w [ci exp[— i J k(x)dx] + c 2 exp[i J k(x)dx]] (10) 

where k — e(x)u>/c, c\ and c 2 are arbitrary constants, and we have assumed that the change in e(x) over one 
wavelength [2ir/k) is sufficiently small compared to |e(x)| (WKB approximation). This gives the usual oscillating 
solution of E z (x, t): 

E z (x, t) w [fc(x)] _ 5 [ Cl ex p[_i( / k(x)dx — tot)] + c 2 exp[i( / k(x)dx — cot)]] (11) 

Jo Jo 
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Since the dielectric has a band gap around the frequency u> 7 these oscillating solutions cannot propagate through it. 
One has to look for exponential or tunneling solutions. In the case of the non-relativistic Schrodinger equation such 
solutions are obtained when the function corresponding to e(x), namely, [E — V(x)], becomes negative. This is not 
possible in electrodynamics because e(x) is never negative. However, it is significant that a general tunneling solution 
can still be found, and is given by 



E^(x, t) ss [ci exp[— / n(x)dx + tut] + C2 exp[ / K,{x)dx — bjt\\ ) 

Jo Jo 



(12) 



with n(x) — ui\J 'e(-ix) j 1 c a real, positive function |13| . This is clearly a solution of the wave equation 

dlEi - *^fl?2S< = (13) 



which is Lorentz invariant as long as e(— ix) is a real, positive Lorentz scalar function. That is guaranteed if e(x, t) is 
a real, positive definite function of the Lorentz invariant variable {x 2 — c 2 t 2 ) in an arbitrary inertial frame. We will 
therefore restrict our discussions to such cases only. 

Notice that the tunneling solution ( |l2| ) is a mapping of the oscillating solution (|ll]) by 

x — > — ix, t — > —it (14) 

Maxwell's equations in vacuo are invariant under this mapping. Maxwell's equations in an inhomogeneous dielectric 
[equations (1) - (4)] are also invariant provided e(— ix) = e(x). But that is certainly not the most general case. 
Assuming that e(x) is an analytic function, one can express it as a Taylor series around x = 0: 

e(x) = e + a nX n (15) 

n 

with the sum positive definite (E). Thus e(—ix) will be complex in general. But, since Im-^/ e(— ix) will give rise to 
oscillating terms in (|l2|), and since the material is assumed to have a band gap around ui, it must vanish. Maxwell's 
equations then get mapped on to equations, such as equation (|l3|), that are still Lorentz invariant and therefore 
acceptable. It is clear from equation (O) that the propagation will be superluminal provided e(— ix) < eo(= 1). This 



is possible, for example, if the dielectric function e(—ix) = (1 + ^2 n a n x n ) < 1 with n such that lxaye{—ixj = and 
E„ a n x n < 0. 

An immediate consequence of the mapping ( |l4| ) is that time-like intervals are mapped on to space-like intervals 
(c 2 t 2 — x 2 ) — * (x 2 — c 2 t 2 ) . Consequently, if e(— ix) < 1, all causally related events get connected by superluminal 
signals. Conversely, it is straightforward to see that superluminal signals (v > c) imply the mapping (|l4| ) , because 



= (x — vt) j yl — v 2 / c 2 = —i(x — vt) j y/ v 2 jc 2 — 1 



t' = {t- vx/c 2 )/y/l-v 2 /c 2 = -iff - vx/c 2 )/y / v 2 /c 2 - 1 (16) 

This is remarkable and important for the interpretation of the experiments showing superluminal tunneling — they 
do not contradict Lorentz invariance and causality. 

It is instructive to look at the difference between superluminal optical tunneling and tunneling of massive particles. 
While tunneling, the energy and momentum of massive relativistic particles are imaginary, as one can easily verify by 
applying the energy and momentum operators on their wavefunction. Thus, the relativistic relation E 2 = p 2 c 2 -\-m 2 ) c A 
gets mapped on to E 2 — p 2 c 2 — m§c 4 , implying tachyons. This does not happen for massless bosons. Nevertheless, 
as we have seen above, tunneling solutions in electrodynamics are also superluminal. 

It is often asserted that according to the special principle of relativity the maximum velocity that a physical signal 
can have is the velocity of light c in vacuum. If that is correct, then the special relativity principle would rule out the 
possibility of dielectric materials of the kind discussed above. That would imply that somehow only dielectrics with 
the property e(— ix) = e{x) can exist physically. Whereas that is not impossible, we find it hard to believe that such 
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a demonstration can indeed be given. On the other hand, if one restricts oneself to the assumptions actually made by 
Einstein, namely the postulate of relativity of uniform motion coupled with the postulate that the velocity of light is 
independent of the motion of the light source, one need only insist on Lorentz invariance as a necessary condition for 
a physical law 14 1. That would leave open the possibility of dielectrics of the kind that would make superluminal yet 
causal signals possible in tunneling modes. 

Interestingly, the dielectrics chosen in the tunneling experiments Q all had variable layers of dielectrics and were 
practically dispersion free. Now, it is well-known that causality and dispersion relations are intimately related |jll| . 
It follows from these dispersion relations that the real part of the refractive index n must vanish for a purely non- 
dispersive material. Hence the velocity of propagation c/n of light through such a material has no upper limit. The 
problem is to produce such materials. The trick is to prepare a medium in such a way that it is inhomogeneous with 
alternative thin layers of high and low refractive indices rii that are all greater than unity (rii > 1) so that it acquires 
a photonic band gap. Then the evanescent wave sees a refractive index < 1, as we have seen, and so propagates 
superluminally without changing shape. 



III. QUANTUM MECHANICAL FORMULATION OF OPTICAL TUNNELING 

We will now show how to give a purely quantum mechanical formulation of this superluminal tunneling behaviour. 
For this we need to use a consistent quantum mechanical formulation of massless electrodynamics using the Kemmer- 
Harish- Chandra formalism outlined in the Appendices. It is clear from this formalism that the classical Maxwell 
fields are components of a ten-component quantum mechanical wavefunction with constraints that reduce the degrees 
of freedom to two. For the tunneling problem, the number of degrees of freedom is further reduced to one, as we 
have already seen. Let the incident finite duration signal be represented by the electric fields (components of the ten 
dimensional unnormalizcd photon wavefunction jip, vide Appendix A) 



E\ — J dkA(k) cos {kx - wt - 4>) - 

J dkA(k) cos (kx + cut + <p) for x < (17) 



1 f x 
E d z = 0(t) Ccxp [- / n(x)dx + (J t\ for < x < d 

\/k(x) ' Jo 



(18) 



E( = 9{t - t)VT J dkA(k) cos [k(x -d)- u(t -t)+x] for x > d (19) 



where A(k) = (l/V2wa 2 ) exp [-(k - k a ) 2 /2a 2 } is real and A(k)dk = 1, J^° oo kA(k)dk = k . R and T are the 
reflection and transmission coefficients, k — u>/c, k(x) = k^^J t(~ix), t is the tunneling or dwell time and 9(t) is the 
step function. (Note that there is no term representing a reflected wavefunction within the tunneling region because 
we are not considering a steady state situation or times t > r.) Accordingly, the dielectric medium is at rest (in 
the sense of being free of any disturbance) before t — and there is no emerging signal at x = d before t = r. By 
matching the wavefunctions smoothly at the boundary x = 0, t = 0, we get 



C = VMO) (1 - VR) cos (20) 
tan</>=4^ = 1 ( 21 ) 



Hence 



7 .i QM VW)(1-^R) cos0 



El = 9(t) v v ' -exp[- / n(x)dx + w t] (22) 

y/K(x) Jo 

The magnetic field in the tunneling region is determined by the analog of (|7|) for the tunneling case and is given by 
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H d = 9(t) — d x E d (23) 

Wq 



Therefore we have (in the WKB approximation) 

c f x 

H d = -e{t)y^^(^-VR)^/KJx)cos(j)—exp[- n(x)dx + w t] (24) 

^0 Jo 

Matching the wavefunctions at the other boundary x — d, t = r gives 

VT = v v a 1 sec x exp _ / n(x)dx + loot] 

Further, matching the derivatives of the wavefunctions at this boundary, one has 



(25) 



tan X =^ ( 26 ) 

/Co 



The velocity operator in this formalism is the 10 x 10 matrix v(3 x = (c/ye(- ix)) 
(Poflx — PxPo)- Thus the Poynting vector can now be calculated, and is given by (see Appendix A) 

S d x = moc^ T 1 p xl 4> = -cE d H d 

n2 



c 2 f x 

0(t)/c(O)(l - VR) 2 cos 2 (f> — exp[-2( / K {x)dx - w t)] (27) 

2^o Jo 



The energy density is given by (see Appendix A) 

£ d = \Vli> = \U-ix)Ef + Hf] 



c 2 f x 

9(t)K(0){l-VR) 2 cos 2 (j)—^K(x)exp[-2( n{x)dx - w t)] (28) 

2w Jo 



One can therefore calculate the velocity of energy transport 



x £ d ^Tixj 

It follows from this that the tunneling time is given by 

dx 



(29) 



(30) 



_ fd. 
which implies 

/ n(x)dx — lu t = (31) 
Jo 



In a hypothetical model in which y/e(—ix) = 1 



r= d -- a -f (32) 
c 3c 



which is always less than the time for passage through vacuum. This supcrluminal effect will be further accentuated 
if one includes higher order terms in x in the expansion of \Jc(— ix) because of the condition J2 n a n x n < stated 
above. 

If one uses the de Broglie-Bohm guidance condition v d = dx/dt, one again obtains the same result for r. These 
results confirm that the energy and so the physical signal indeed propagates supcrluminally while tunneling. 
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IV. CONCLUSIONS 



In conclusion we would like to emphasize precisely the significant new result that we have obtained. Since there 
has been much discussion and some controversy in the literature regarding superluminal effects and their causality, 
let us summarize the situation as we see it. 

The materials used for observing superluminal effects have been generally termed "ultrarefractive" jlq| . Near the 
edges of a transmission gap the effective permitivity can become close to zero. Consequently, surprising effects can be 
observed on light transmitted and reflected by such materials, such as superluminal velocities as well as enlargement 
and splitting of the transmitted beam. 

In one type of process the effects are results of anomalous dispersion, i.e., anomalous variation of the permitivity 
with wavelength. Although the 1914 analysis of Sommcrfeld and Brillouin clearly established that superluminality 
in such cases cannot be Einstein causal and is only apparent, it has recently been argued that this need not be the 
case for physical signals that are of finite duration and extent because a responsive signal travelling backward in time 
in such a case cannot arrive before the primary signal is generated, thus preserving the causality principle Q. Our 
paper does not deal with this type of phenomena. 

The second type of process involves tunneling in one (or two) dimensions through a narrow band gap, and it is only 
this type of phenomena (ID tunneling) that we have addressed. The theoretical discussions of such phenomena have 
so far been based purely on an analogy between the non-relativistic Schrodinger equation and the Helmholtz equation 
leading to an effective refractive index n(x, y, z) — {2m[E — V(x, y, z)\\ 1 ^' 2 c/fujj which is imaginary in any region where 
E < V [^,^). This mechanism is, in reality, not applicable to photons, as we have mentioned earlier and as Chiao and 
Steinberg admit in their review article fl6|| . To take a definite stand on an issue such as superluminal propagation and 
causality, analogies are not reliable in our opinion, and one must use a proper theory, namely a consistent relativistic 
quantum mechanical formalism for photons Q . We have used this formalism to carry out explicit calculations for the 
tunneling of a finite width photon wave-packet incident normally on a ID photonic barrier. (Note that in this sense also 
our result is new because total internal reflection in optics occurs only for non-zero critical angles of incidence.) Our 
analysis clearly shows that genuine Einstein causal superluminal propagation can occur only if the tunneling medium 
is inhomogeneous on the scale of the wavelength and Ime(-ix) — 0. This follows simply and very generally from the 
fact that points on the light cone remain on the light cone under the mapping ( |T4| ) which takes propagating solutions 
to tunneling solutions. Therefore, the only way to get genuine superluminal signals is to have an inhomogeneous 
dielectric function e(x) > 1 that is mapped to e(— ix) < 1 with Ime(-ia;) = to ensure Lorentz invariance of the wave 
equation ([l^). This argument obviously holds for both classical and quantum light, and is consistent with dispersion 
relations and causality JTT[ ] . 

Such materials have been used in actual experiments |^6| , ^7| . They involve tunneling at near normal incidence 
through band gaps excited in periodic dielectric structures. These band gaps arise from Bragg reflections from the 
periodic structure, leading to an evanescent decay of the wave amplitude when the frequency is within the forbidden 
band gap at the first Brillouin zone. It should be noted that such periodic structures are non-dispersive so that 
the tunneling wave-packets that are tuned to midgap remain essentially undistorted upon transmission through the 
barrier, though much attenuated in amplitude flq| . 
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VI. APPENDIX A 

Until recently, no consistent quantum mechanical formalism existed for relativistic bosons below the threshold for 
pair production and annihilation. Relativistic quantum mechanics can only be consistently formulated provided there 
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exists a conserved four-vector current whose time component, to be identified with the probability density, is positive 
definite. Unfortunately, the conserved charge vector current for relativistic spin and spin 1 bosons does not have this 
property. Moreover, the charge current vanishes for neutral particles like the photon. However, it has now been shown 
H that a conserved four-vector current with a positive definite time component does exist for relativistic bosons, and 
is associated, not with the charge current but, with the flow of energy. This formulation is based on the first-order 
Kemmer equation |^| 

(ihP^d^ + m c) tp = (33) 

where the matrices (3 satisfy the algebra 

/9 M A, P\ + P\ Pv = P» 9v\ + Px g»n ■ (34) 

The 5x5 dimensional representation of these matrices describes spin bosons and the 10 x 10 dimensional represen- 
tation describes spin 1 bosons. Multiplying ([33]) by Po, one obtains the Schrodinger form of the equation 

i h — — = [—ihc (3idi — to c 2 (3q } ip (35) 
at 

where Pi = P Pi — Pi P . Multiplying ( |33| ) by 1 — /3q, one obtains the first class constraint 

ihf3 i p 2 d l 4> = -m c(l- [3l)ip. (36) 
It implies the conditions divD = — (mQc/h)Ao and B = curLA if one takes 



V> T = (l/V m c 2 )(-D x ,-D y , -D z , B x , B y ,B z ,-moA x ,-moAy,-moA z ,mA )) (37) 

The reader is referred to Ref. || for further discussions regarding the significance of this constraint. 

If one multiplies equation (|35| ) by ip^ from the left, its hermitian conjugate by if> from the right and adds the 
resultant equations, one obtains the continuity equation 

d W+MiM = o. (38) 



ot 



This can be written in the form 



where 



0" ©mo = (39) 



®»u = -~m Q c 2 ^(f3^ + pvpft - g^)if} (40) 
(with ip = ip'r)Q, t]q = 2/3q — 1,T)q = 1) is the symmetric energy-momentum tensor, and 

e 00 = -toqc 2 ^ V < o (4i) 



Thus, it is possible to define a wave function (f> = y/moc 2 / E ip (with E — — J O o dV ) such that <ft <f> is non- 
negative and normalized and can be interpreted as a probability density. The conserved probability current density 

is SjU = -q^o/e = (tfcf), -ft pi <p). 

Notice that according to the equation of motion (j35"[), the velocity operator for massive bosons is c/3j. 
The theory of massless spin and spin 1 bosons cannot be obtained simply by taking the limit too going to zero 
because of the l/y'mo factor in ip. One has to start with the equation 1 1 1 

ihP^ip + mocjii = (42) 

where 7 is a matrix that satisfies the following conditions: 
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7 2 - 7 (43) 
lP»+M = Pi>- ( 44 ) 
This equation can be derived from the gauge invariant Lagrangian density 

c = - j [Pfaw - W^d^} + ^f^hl> (45) 
Multiplying ( ff2| ) from the left by 1 — 7, one obtains 

/3 M a"( 7 ^) = 0. (46) 
Multiplying (^) from the left by d\ f3 x fl" , one also obtains 

a A /?A/3,(7^) = ^( 7 ^). (47) 

It follows from © and © that 

□ (7V) = (48) 

which shows that 7 ^ describes massless bosons. 
The Schrodinger form of the equation 

in^^- = -ihckdi{-iijj) (49) 

and the associated first class constraint 

ihfcPl flj^ + Tn o c(l-/3g)7V' = (50) 

follow by multiplying ( f42| ) by /3o and 1 — /3q respectively. Equation (|4^) implies the Maxwell equations curLE = 
-(jj,/c)dtH and curLff = (e/c)d t E if 



7^ T = (l/y/moc?)(-D x , -D y , -D 2 , B x , B y , B z , 0, 0, 0, 0) (51) 
The constraint (|50| ) implies the relations divi? = and B — curLA. The symmetrical energy-momentum tensor is 



and so the energy density 



2 

&»» = - 1 ^$Wv+PuPv-9 v »)'Vl> (52) 



£ = -Soo = ^ V> f 7^ = \ [E.E + B.B] (53) 



is positive definite. The rest of the arguments are analogous to the massive case. 
The Bohmian 3-velocity Wj for massless bosons can be defined by 



(54) 



Notice that in relativistic quantum mechanics the Bohmian velocity is not defined through the gradient of the phase 
as in non- relativistic quantum mechanics but in terms of the energy flux current. 

Neutral massless vector bosons are very special in quantum mechanics. Their wave function is real, and so their 
charge current = ip T (3^ ji/j vanishes. However, their probability current density does not vanish. Furthermore, 
the Poynting vector turns out to be given by 

Si = tooc 3 ^/^ = 4E x H]i (55) 

One might wonder about the significance of the mass parameter mo for massless electrodynamics. It is necessary 
for a consistent quantum mechanical formalism for dimensional reasons and drops out of all physical results because 
of the operator 7. It can be altogether eliminated in favour of the intrinsic parameters in the theory, namely c, H, the 
frequency lj and the spin multiplicity s. 

The representations of the Kemmer-Duffin-Petiau (3 matrices used in this paper are given in Appendix B. 
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VII. APPENDIX B 
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